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1 Introduction 

The majority of special functions and orthogonal polynomials introduced during the last decade are 
associated with Lie groups or their generalizations. In particular, special functions of mathematical 
physics are in fact matrix elements of representations of Lie groups [49] and recent multivariate 
generalizations of classical hypergeometric orthogonal polynomials are based on root systems of 
simple Lie groups/algebras [3, 10, 20, 30, 31, 32, 41, 47]. In this connection a number of elegant 
results in theory these families of orthogonal polynomials, such as explicit (determinantal) compu- 
tation of polynomials [22, 23, 46] and Pieri formulas [25, 48], were obtained, see also [24, 26, 21] 
and references therein. 

The main purpose of this article is to construct orthogonal polynomials in n variables based on 
orbit functions related to simple Lie groups of rank n (see Section 2.2 for definitions). As far as we 
deal with the functions invariant /skew-invariant under the action of the corresponding Weyl group 
the obtained polynomials appear as building blocks in all multivariate polynomials associated with 
root systems. Unlike Gram-Schmidt type orthogonalization of the monomial basis with respect to 
Haar measure [10, 32, 41] or determinantal construction of polynomials [22, 23, 46] we make profit 
from decomposition of products of Weyl group orbits and from basic properties of the characters 
of irreducible finite dimensional representations. 

Relating the polynomials to the Lie groups, by means of either the characters or closely related 
to it the Weyl group orbit functions (C- and 5- functions) , allows one to carry over powerful results 
of the accomplished theory of the compact simple Lie groups as properties of the polynomials. Let 
us point out the following properties: 



(a) The domain of orthogonality F of characters and/or orbit functions coincide in the real 
Euclidean space M". The domain F is known [1] for all G. The polynomial substitution of 
variables transforms F into the domain of orthogonality F for the polynomials. Examples of 
rather bizarre shape of F of the Lie groups C2 and G2 are shown in [12] and [38] respectively. 

(b) Orthogonality of characters, C-, and S'-functions, when integrated over F, is known [37]. 
Prom it follows directly the orthogonality of the polynomials when integrated over F. 
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(c) Discretization and discrete orthogonality of characters and C-functions of [36] extends directly 
to S'-functions [37] . For any G it is done on the fragment of a lattice in F of any chosen density. 
Lattice points in F become discrete set of points in F after the polynomial substitution of 
variables. 

(d) Congruence classes of characters, orbit functions and of the polynomials are a practically 
useful in computing with these objects [29]. They arise from the action of the center of the 
corresponding compact simple Lie group. 

Other useful features of orbit functions [15, 16] should be reflected in properties of polynomials 
originating from simple Lie groups: in particular, orbit functions are eigenfunctions of differential 
operators, the Laplace operator being one of them. C- and S'-functions are solutions of the Neumann 
and Dirichlet boundary value problem respectively. However, our primary objective at this stage 
is to establish a constructive method for finding orthogonal multivariate polynomials, indeed, for 
actually seeing them. 

Our method is purely algebraic and we propose three different ways to transform a C- or S-orbit 
function into a polynomial. 

The first one substitutes for each multivariablc exponential term in an orbit function a mono- 
mial of as many variables (3.1). In ID this results in Chcbyshcv polynomials written as Laurent 
polynomials with symmetrically placed positive and negative powers of the variable; and in the 
case of A2 our results coincide with those from [19]. 

The second method, the 'truly trigonometric' form, is based on the fact that, for many simple 
Lie algebras (see the list in (3.3) below), each C and S'-orbit function consists of pairs of exponential 
terms that add up to cither cosine or sine. Hence such a function is a sum of trigonometric terms. 
For the Chebyshev polynomials we obtain in this way their trigonometric form. Note from (3.3) 
that this method does no apply to the groups A„ for n > 1. 

This paper focuses on polynomials obtained by the third substitution of variables, mimicking 
Weyl's method for the construction of finite dimensional representations from n fundamental rep- 
resentations (see for example [5], Supplement, Section 5). Thus the C-polynomials have n variables 
that are the C-orbit functions, one for each fundamental weight uij, see (2.1). This approach re- 
sults in a simple recursive construction that allows one to represent any orbit function/monomial 
symmetric function in non Laurent polynomial form. 

In addition to the general approach and associated tools we present a lot of explicit and prac- 
tically useful data and discussions, namely in Appendix A we compare the classical Chebyshev 
polynomials (Dickson polynomials) and orbit functions of Ai with their recursion relations. Not 
quite standard is addition of the ^i-character formula (A. 4), linking C- and ^-polynomials. Suit- 
ably normalized, the Chebyshev polynomials of the first and second kind coincide with the C- and 
^-polynomials. A table of the polynomials of each kind is presented. Appendices C, B, and D 
contain respectively the recursion relations for polynomials of the Lie algebras A2, C2 and G^- In 
Appendix E recursion relations for A3 and generic recursion relations for and C3 polynomials 
of both kinds are listed together with useful tools for solving these recursion relations, namely the 
formulas determining the congruence class of polynomials and the dimensions of the irreducible 
representations . 

2 Preliminaries and conventions 

This section serves to fix notations and terminology and to recall the definitions and some of the 
properties of orbit functions. Additional details can be found for example in [1, 11, 13, 14, 15, 16, 
17, 18]. 
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2.1 Notations 



Let M" be the Euclidean space spanned by the simple roots of a simple Lie group G. The basis 
of the simple roots and the basis of fundamental weights are hereafter referred to as the a-basis 
and (jj-basis respectively. The two bases are linked by the Cartan matrix M of G. In matrix form 
that is 

a = Muj, u; = M-^a, M = (M.^) = (^^^^^ , j,k = {1,2, ... ,n}. 

Bases dual to a- and cj-bases are denoted by a- and cD-bases. In addition one uses {ei, . . . , 6^}, 
the orthonormal basis of M". Note, that {aj,cjk) = {otj,0Jk) = i^jj^k) = Sj^, where (• , •) is the 
inner product and Sjk is the Kronecker delta. 

The root lattice Q and the weight lattice P of G are formed by all integer linear combinations 
of the a-basis and w-basis, 

Q = Zai+Za2 + --- + Zan, P = 'Lui + 1.^2 + ■ ■ ■ + 'LoJn- (2.1) 

In the weight lattice P, we define the cone of dominant weights P~^ and its subset of strictly 
dominant weights P~^~^ 

P Z) P+ = Z^^uji + ■■■ + 'L-^Un D P++ = Z>°a;i + • • • + Z>°w„. 

Hereafter W = W{G) is the Wcyl group, i.e., the finite group generated by reflections in (n— 1)- 
dimensional hyperplanes orthogonal to simple roots, having the origin as their common point, and 
referred to as elementary reflections Vj, j = 1, . . . ,n. The orbit of W containing the (dominant) 
point A G P+ C M" is written as Wx. The size of Wx is denoted by \Wx\, it is the number of points 
in Wx. 

The fundamental region F{G) C M" is the convex hull of the vertices {0, . . . , ^}, where qj, 
j = l,n are comarks of the highest root ^, i.e., ^ = qiai + • • • + QnCtn- 

2.2 Orbit functions and the character 

An orbit function of n variables is the set of distinct points in M" generated by the action of W{G) 
on A. 

Definition 1. The C-function Cx{x) is deflned as 

Cx{x):= e^'''^^'''\ xeW, AgP+. 

Definition 2. The S'-function Sx{x) is deflned as 

Sx{x) := Yl (-l)P('')e2'^'<'^'=^>, xeW, AG P++ . 

where p(/x) is the number of elementary reflections necessary to obtain /x from A. 

The same fi can be obtained by different successions of reflections, but all shortest routes from 
A to ^ will have a length of the same parity, so S'-functions are well deflned. 

In this paper, we always suppose that X, fi ^ P are given in w-basis and x G M" is given in 

n n n 

a-basis, namely A = ^ ^j^j^ A* = f^j^j^ ^jiP'-j ^ ^ 2; = ^ XjOtj, xj G M. Hence the orbit 
j=i j=i j=i 

functions have the following forms 

n 

27Ti Y2 f^j'-^j 

Cx{x)=Y^ -1 = ^ J|e2-'^^-^ (2.2) 

fM&Wx msWa i=i 
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Sxix)= ^ i-lf^f'^e = ^ JJe^"^''^^^. (2.3) 

neWx fi&Wx j=i 

The number of exponential summands that form an irreducible S'-function Sx is always equal to 
the order of the corresponding Weyl group, because the function is non-zero only if the stabilizer 
of A in is the identity element. The number of exponential summands that form an irreducible 
C-function Cx is equal to the order of W divided by the order of the stabilizer of A in W. Note 
that in the 1-dimensional case, C- and S'-functions are respectively cosine and sine functions up to 
a normalization. The special case, C(o,o,...,o) = 1 G, occurs if A is the origin of M". 

Occasionally it is useful to scale up Cx of non-generic A by the stabilizer of A in T^. It is done by 
somewhat modifying the Definition 1. Rather than summation over the elements of Wx, one should 
sum over all the elements of W. Then A would be counted as many times as there are elements 
w & W stabilizing A. That is the elements with the property wX = X. 

There is a fundamental relation between the C and S-orbit functions for simple Lie group G of 
any type and rank, called the Weyl character formula. For the character xx{^) of the representation 
of G with the highest weight A, it is the following: 

XA(x) = %f^ = E"^M^M(^)' ^^^^ A,^GP+. (2.4) 

n 

Here p is the half sum of the positive roots of G, which is known to be for any G given by p = Wfe. 

k=l 

The positive integer is the Kostka number (the multiplicity of the dominant weight ji in the 
weight system of the irreducible representation with the highest weight A). A basic algorithm in 
representation theory [34] is used for computation of multiplicities. The multiplicity of the highest 
term in (2.4) is always one, i.e. m\ = 1 in all cases. A suitable ordering of the weights A G P+ makes 
the matrix (m^) triangular (see Tables in [2]) with all diagonal entries equal to one. Such a matrix 
can be easily inverted, so that any C^{x) can be written as a linear combination of irreducible 
characters with coefficients that are integer but not all positive. 

C^{x) = Y,{ml)-\x{x) = S;^Y.^ml)-^Sx+p, x G M^ A,/xGP+, (2.5) 

\ A 
where (??t.^)~^ are elements of the triangular matrix inverse to (m^). 

Example 1. This example illustrates the relations (2.4) and (2.5) between polynomials C and x 

of the Lie group/Lie algebra of type G2. 

Tables 1 contain values of and {m^)~^ for the five lowest dominant weights of G2. It is thus 
possible to write any of the lowest five characters as the linear combination of C-polynomials (2.4), 
as well as the inverse relation, namely writing C-polynomials in terms of the characters (2.5). 
In particular, using the entries in the last column of the first table, we obtain 

X(i,i) = 4C(o,o) + 4C(o,i) + 2C(i^o) + 2C(o,2) + (2-6) 
Similarly according to (2.5), one reads the last column of the second table as follows, 

C(i,i) = 2x(o,o) - 2x(o,2) + (2-7) 

A useful verification: The equalities in (2.6) and (2.7) must be maintained when x\ and C^ are 
replaced by the dimensions of the representation (see (D.l)) and by the size of the W-orbits (see 
Table 1) respectively. 
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^) ^ for the lowest 5 dominant weights of G2. The first 



Table 1: Tables of multiplicities and 
row contains A, the first column of each table contains The last column contains the size of the 
Weyl group orbit of the dominant weight /x. The last row contains the dimensions of the irreducible 
representations of G2 with the highest weight A. 



2.3 Properties of orbit functions 

The rank of the underlying semisimple Lie group/algebra is the number of variables of the three 
families of orbit functions. In general, C- and S-functions are finite sums of exponential functions, 
therefore they are continuous and have continuous derivatives of all orders in M". 5- functions are 
antisymmetric with respect to the (n — l)-dimensional boundary of F. Hence they are zero on 
the boundary of F. C-functions are symmetric with respect to the (n — l)-dimensional boundary 
of F. Their normal derivative at the boundary is equal to zero (because the normal derivative of 
a C-function is an ^'-function). Symmetry and other properties of the orbit functions are reviewed 
in [15, 16, 17]. 

Through the use of a method of orbit functions it is possible to solve eigenvalue problems on 
the fundamental domain for any compact simple Lie group. Namely it was shown [15, 16] that 
C\{x)- and S'A(x)-orbit functions are eigenfunctions of the n-dimensional Laplace operator on the 
simplexes, which are fundamental domains of compact simple Lie groups, with the Neumann or 
Dirichlet boundary value conditions. The Laplace operator has the same eigenvalues for every 
exponential function summand of an orbit function, the eigenvalue equals to — 47r(A, A). Indeed, all 
points of the orbit containing A are equidistant from the origin. 

For any complex square integrable functions (f){x) and tlj{x), we define a continuous scalar product 

(0(x),V(x)) := / (/.(x)^dx. (2.8) 

JF 

Here, integration is carried out with respect to the Euclidean measure, the bar means complex 
conjugation, and x € F, where F is the fundamental region of W. 

Any pair of orbit functions from the same family is orthogonal with respect to the scalar prod- 
uct (2.8) on the corresponding fundamental region [11], namely 

{Cx{x),Cx'{x)) = \Wx\ ■ \F\ ■ Sxy, (2.9) 
{Sx{x),Sy{x)) = \W\-\F\-6xy, (2.10) 

where Sw is the Kronecker delta, |Ty| is the order of Weyl group, \W\\ is the size of Weyl group 
orbit, and \F\ is the volume of fundamental regions (formulas for \F\ are found in [11]). 

The functions of the families C and S are complete on the fundamental domain (completeness 
follows from the completeness of the exponential functions). 
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3 Multivariate orthogonal polynomials corresponding to orbit functions 



In this section wc describe three substitutions of variables 

that transform the C(x)- and >S'(a;)-orbit functions into their polynomial forms in n variables. 

It directly follows from the orthogonality of the orbit functions (2.9) and (2.10) that such polyno- 
mials are orthogonal on the domain F, where F is the image of the fundamental region F under the 
transformation T. The corresponding integration weight function can be found using the Jacobian 
det~^ (^^iy) transformation T. 

3.1 Polynomials obtained by exponential substitution 

The first substitution of variables is rather straightforward 

Xj := e^'^*^^' , Xj eR, j = 1,2,..., n. (3.1) 

n 

Polynomial summands are products Yl -^f^ > where G Z are components of the orbit points 

j=i 

relative to a suitable basis. Under this transformation orbit function, Cx{x) and Sx{x), given 
by (2.2) and (2.3), become Laurent polynomials in n variables Xj, where j = {1, 2, . . . ,n}. 

The exponential substitution polynomials are complex-valued in general, admit negative powers, 
and have all their coefficients equal to one in C-polynomials, and 1 or —1 in S'-polynomials. 

Common trigonometric identities can be viewed as identities between C- and S-orbit functions 
of one variable. It is likely that identities between C- and S-orbit functions of more than one 
variable could also be found. 



3.2 Polynomials obtained by trigonometric substitution 

The M^-orbits of many simple Lie groups (but not all!) have an additional property that admits 
a truly trigonometric substitution of variables. In order that Wx{L) is such an orbit we have to 
have the following, 

±HeWx{L) for all fjLeWx{L). (3.2) 
The pair of corresponding terms of the function of Wx{L) can be combined: 

^2ni{^,x) ^ ^-2ni{^,x) ^ 2cos(27r(^,.x)) e Cx{x), 

_ ^ 2isin(27r(Ai,.x)) € Sx{x), 

so that Cx{x) and Sx{x) become linear combinations of cosines and sines. Using common trigono- 
metric identities, these cosines and sines can be expressed through the lowest ones, which are thus 
chosen as the new polynomial variables. 

Obviously the important question is when (3.2) is valid. In Lie theory, the answer is known for 
representations, and extends without reservation to I^-orbits. The property that assures validity 
of (3.2) is self-contragrediency of the representation of the highest weight A (see for example [33]). 

All VF-orbits of the following Lie groups have the property (3.2), hence admit the trigonometric 
substitutions: 

Ai, Bn (n > 3), Cn (n > 2), Dan {n > 2), Eg, F4, G2. (3.3) 

Note that in the case of Ai, this is precisely the trigonometric substitution made for Chebyshev 
polynomials of the first and second kind. 

The representations/orbits of the remaining Lie groups that also have the property (3.2) are 
listed in [33]. 
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3.3 Polynomials obtained by recursive method 

In the previous subsection, the substitution of variables transforming orbit functions into polynomi- 
als can be seen as generalization of the trigonometric variables of classical Chebyshev polynomials. 
Here, the substitution of variables is also a generalization of the trigonometric variables of classical 
Chebyshev polynomials, but of a different kind: It works uniformly for simple Lie algebras of all 
types, not only those of (3.3). 

The fundamental weights of a simple Lie algebra L„ of rank n are the basis vectors u^, k = 
1, 2, . . . , n, of the w-basis of W^. We choose C-functions of the n fundamental weights as the n new 
variables: 

Xj:=Xj{x):=C^.{x), j = l,2,...,n, x G R"; (3.4) 

completed by one more variable, the lowest S'-function, 

n 

S:=Sp{x), xGM", p=il,l,...,l) = ^Um, (3.5) 

m=l 

which cannot be constructed by multiplying the other variables. 

The recursive construction of C -polynomials begins by multiplying the variables Xj and C- 
functions and decomposing their products into sums of C-polynomials. A judicious choice of the 
sequence of products allows one to find ever higher degree C-polynomials. 

First, generic recursion relations are found as the decomposition of products XjC(^ax,a2,...,an) 
with 'sufficiently large' ai,a2, . . . ,an (i.e. all C-functions in the decomposition should correspond 
to generic points). Then the rest of necessary recursions additionaV) are constructed. An efficient 
way to find the decompositions is to work with products of Weyl group orbits, rather than with orbit 
functions. Their decomposition has been studied, and many examples have been described in [7]. 
Note that these recursion relations arc always linear and the corresponding matrix is triangular. 
The procedure is exemplified below for L = A2, A^, C2, and G2. Also generic recursion 
relations are shown for and C3. 

Results of the recursive procedures can be summarized as follows. 

Proposition 1. Any irreducible C -function and any character x\ of a simple Lie group G can be 
represented as a polynomial of C -functions of the fundamental weights ui, ... ,ujn, i.e. a polynomial 
in the variables Xi,X2, . . ■ , X^. 

Proof. Let us fix the dominance order (natural partial order) on the weight lattice: 

A ^ /Lt (A - /x) G Q+, 

where = ^ Na, A+-positive roots. 

aeA+ 

Then we say that 

As soon as the above ordering is fixed we can order our polynomials Ca(X), A = Aia;i+- • --t-A^w^. 
Then the above proposition follows from the method of induction with respect to weight lattice 
point A. 

The basis of induction is formed by the additional recursions, i.e., by all such expansions of 
the products XjCx that contain at least one orbit function C^, where A have at least one zero 
coordinate in w-basis. 
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As fax as A is given in w-basis and x is given in d-basis the expansion of the product Xj{x)C\{x) 
is equivalent to the expansion of the Weyl orbits product W^.W\. Here exists unique A' € W\ such 
that V/Lt € Wx we have /j, ^ X'. Then the general form of the recurrence relation is following 



XjCx = C(y^,A^,...A;.+i,...x) + ^(C^A'), i = 1, 2, . . . , n. (3.6) 

Where L(C^a') is Z-linear combination of orbit functions C^, with /x =<; A'. The formula (3.6) 
implies the induction step for each j = 1,2, ... ,n. 

Therefore we proved that C\ is the polynomial of C-functions of the fundamental weights 
uji, . . . ,u}n- The same statement for xx follows from the fact that character is the linear com- 
bination of certain C-functions. □ 



The recursive construction of S-polynomials starts by multiplying the variables S and Xj and 
decomposing their products into sums of S'-polynomials. Again a judicious choice of the sequence 
of products allows one to find ever higher degree S'-polynomials. The procedure is exemplified 
below for L = A2, C2, and G2. However, let us point out that, the higher the rank of the underly- 
ing Lie algebra, the recursive procedure for S'-polynomials becomes more laborious, in comparison 
with the similar procedure for C-polynomials. This is caused by the presence of negative terms 
in S'-polynomials, and by the fact that all S-functions have the maximal number of exponential 
summands. In calculating the decompositions of products of polynomials there are frequent can- 
celations of terms. 

Fortunately, there is an alternative to the recursive procedure for S-polynomials. Once the C- 
polynomials have been calculated, they can be used in (2.4) for finding S'-polynomials as sums of 
C-polynomials multiplied by the variable S. 

Remark 1. It is important to note that in practice, polynomials should be used instead of S^. 
This is caused by the fact that the polynomial Sx depends on n + 1 variables S, Xi, . . . , Xn but 
it is defined for n-dimensional Euclidean space. The fraction allows us to avoid this confusion. 
The function S(X) has no zeros inside F and at the boundary of F both S and Sx have zero 
simultaneously, so the fraction ^ is well defined. In fact this means that we are working with the 
character instead of S'-orbit function. 

Remark 2. There are two easy and practical checks on recursion relations applicable to all simple 
Lie algebras. The first one is the equality of numbers of exponential terms in S- or C-functions on 
both sides of a recursion relation (the numbers of exponential terms are calculated using the sizes 
of Weyl group orbits) . 

The second check is the equality of congruence numbers The character xx^ tlie representation 
A, and the Cx- and S'^-polynomials, as well as all the points in one Weyl group orbit, can be assigned 
to a congruence class specified by the congruence number #(A), which is the number (label) of the 
coset containing A in the lattice P with respect to the root lattice. In other words, our congruence 
numbers # are the elements of the quotient group P/ Q. The number of possible distinct values of # 
is equal to the order of the center of G. The congruence numbers add up during multiplication. In 
particular, all the C-polynomials in (2.4) must be from the same congruence class. More generally, 
all the terms in the decomposition of a product of the polynomials belong to the same class. This 
criterion is trivial for three of the simple Lie groups: C2, F4 and Eg, since they have only one 
congruence class. 

The number of exponential terms in a character xxix) is equal to the dimension of the represen- 
tation with the highest weight A. Hence the number of terms in each of the C-polynomials in (2.4) 
have to add up to the dimension of the representation. 
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Table 2: Number of points in VF-orbits for the specific Lie algebras considered in this paper. The 
orbits are labeled by their dominant weights A, where the symbol stands for any positive integer. 

4 Relations between orbit function polynomials and Macdonald polynomials 

In this section we establish connections between orbit functions and known multivariate orthogonal 
polynomials. In particular it will be shown that C-polynomials play the role of building blocks in 
the theory of symmetric and antisymmetric polynomials (Schur polynomials, Jacobi polynomials of 
many variables, Macdonald symmetric polynomials). 

4.1 Monomial symmetric polynomials 

Orbit functions C\ are a certain modification of monomial symmetric polynomials 

mx{y) = E = E yi'y2" ■■■vn", a g p+. (4.i) 

This relation follows from the substitution (3.1) introduced in the previous section (here yj = Xj, 
i = l,2,...,n). 

For studying symmetric orthogonal polynomials one usually replaces y^ by {e^ is considered 
as a function on Euclidean space e^(x) = e^^'^^ = e'^^^^"' i-^*"^"). Then we obtain modified C-orbit 
functions also called monomial symmetric functions 

mx{x) = E e^'''''^ = E e'^i^i+-+'^"^". (4.2) 

The functions mA(2;) are also eigenfunctions of the Laplace operator. Note that if wc take integral 
orthogonal coordinates mi, 7712, . . . , m„ in the An case in such a way that mi > m,2 > • • • > m„ > 0, 
then Laurent polynomials m\{y) and mA(a:;) turn into usual (not Laurent) polynomials. 

Jacobi polynomials in one variable are well-known orthogonal polynomials of the theory of 
special functions and multivariate Jacobi polynomials are symmetric Laurent polynomials, which 
are defined by means of ttia, A G P+, see [9, 8]. 

4.2 Macdonald symmetric polynomials 

Wc consider Macdonald symmetric (Laurent) polynomials, which are a quantum analogue of Jacobi 
polynomialsa and constructed by means of monomial symmetric polynomials. 

Let A denote the group algebra over M of the free Abelian group P, then Weyl group W acts 
on P and A. Further let A^ denote the subalgebra of VF-invariant elements of A and monomial 
symmetric functions rflA defined in previous section form a basis of A^. 

We introduce a real number q, ^ q < 1 and with every root a from the root system R oi G 
we associate a variable tee such that t^ = twa-, w eW. 
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Let Cq^tc be the field of rational functions in q and ta and let Cq^ta [^^^ i • • • > e^"] denote the set of 
Laurent polynomials in functions from the previous section e^^ , e^^ , ■ ■ • , e^" with coefficients from 
C{q,ta). 

The constant term ao of a function J2 '^ac'^ we denote by aAe^jo- 

A A 

For f,g& ^q,ta [c^^ ) • • • ) c^"] we determine an inner product 

oo -. 2i a 

(/,5W = n-M/5A]o, where A = 11 i _7 J,. ^ 

and the bar over g denotes the linear involution determined by = e~^. 

The Macdonald polynomials Px, X E P'^ are uniquely defined by the following theorem proved 
by L Macdonald [32]. 

Theorem 1. There exists a unique family P\ G Cq^t^[e^^ , e^"^ , . . . , e^"]^ , A G P+, satisfying the 
conditions 

1) Px=mx + Y.«^ «AeC,,t„, 

At<A 

2) (PA,^^),,t„ =0, if 

In other words the Macdonald polynomials arc obtained via orthogonalization of the basis for 
. The orthogonality property of the Macdonald polynomials is proved by showing that the 
Macdonald polynomials are eigenvectors for an algebra of commuting self adjoint operators with 
one-dimensional eigenspaces, and using the fact that eigenspaces for different eigenvalues must be 
orthogonal. 

Replacing e^^ by yj, j = 1, 2, . . . , n, in Px we obtain (for each fixed values of q and ta) orthogonal 
symmetric polynomials which are called Macdonald symmetric polynomials. 

Replacing yj by e^"^'^^ , j = 1, 2, . . . , n, in Macdonald polynomials we obtain orthogonal functions 
which arc finite linear combinations of C-orbit functions. 

Some special values of q and ta reduce Macdonald polynomials to C- and S- orbit functions: 

• If = 1, then independently of q we have Px{y) = ^xiv) = Cx{x)\y=2mx- 

• UyaeR ta = q, then Px{y) = xxiv) = ^\y=^- 

Corollary 1. Polynomial forms of C- an S-functions introduced in Section 3 by substitution (3.1) 
are partial cases of the Macdonald symmetric polynomials. Polynomial forms of orbit functions 
obtained by the trigonometric substitution and substitution (3.4) are equivalent to these Macdonald 
symmetric polynomials. 

Theorem 1 and Proposition 1 imply the following conclusion. 

Corollary 2. Macdonald symmetric polynomials can be represented as polynomials in fundamental 
C -functions: Xi = C^^ , • • • , Xn = Cu,„ ■ Such an approach produce not Laurent polynomials, but 

polynomials with all nonnegative powers. 

Remark 3. All C- and S'-orthogonal polynomials described in Section 3 inherit from orbit func- 
tions important discretization properties. An uniform discretization of these polynomials follows 
from their invariance with respect to the affine Weyl group of G and from the well-established 
discretization of the fundamental region F{G), see [11] for details. 

It is worth to mention that behind each polynomial we have special functions (orbit functions). 
It makes our orthogonal polynomials richer and gives a number of advantages, one of them is 
cubature formula introduced in [35]. 
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Due to the relation between orbit functions and Macdonald polynomials the mentioned dis- 
cretization method can be adopted and applied to Macdonald polynomials. 

The orbit functions do not require any sophisticated theory to work with. In general, for a given 
group G, most of the properties of orbit functions come from the Weyl group and the irreducible 
character of G or from the basic notions of the representation theory. This fact makes multivariate 
polynomials obtained from orbit functions attractive and handy for a wide range of applications, 
such as discrete and interpolation problems, cubature formulas [35], reduction of polynomials [40], 
etc. 

For the application reason below we present recurrence relations and lowest polynomials for the 
simple Lie groups Ai, A2, A^, C2, G2, C3, B3. All, but two last cases, contain both generic and 
additional recursions. For groups C3 and B3 we obtain only generic recurrence relations and lowest 
polynomials necessary to solve them are available in [40]. 

The content of the following sections is also motivated by the fact that calculation of additional 
recurrences is not suitable for complete computer automatization. However, as soon as additional 
recurrences were obtained, all other calculations concerning polynomials and their applications 
become very algorithmic and easily can be done by computer algebra packages for Lie theory. In 
following sections it is also shown in all details how the recursive algorithm proposed in Section 3 
works. 

5 Discussion 

• Recent years Lie groups became a backbone of a segment of the theory of orthogonal poly- 
nomials of many variables and monomial symmetric functions/orbit functions became build- 
ing blocks of multivariate orthogonal polynomials. In particular, this idea is illustrated by 
orthogonal polynomials of orbit functions, that are the special cases of the Macdonald poly- 
nomials [32], which in turn are the special cases of more general notion of Koornwinder 
polynomials [20]. 

• There is an alternative way to our construction of the polynomials in all but in the An cases. 
The crucial substitution (3.4) can be replaced by 

Xk ■■= Xojki^), A; = 1,2, ... ,n. (5.1) 

In (5.1) the variables are characters of irreducible representations with highest weights given 
as the fundamental weights, while in (3.4) the variables are C-functions of the fundamental 
weights. Only for An the two coincide, Ca;^. (x) = Xuj^ (•^) all A; = 1, . . . ,n and for all x S M". 
Already for the rank two cases other than A2 there is a difference. Indeed, (5.1) reads as 
follows, 

C2 : Xi = XuiA^) = Cuiix), X2 = Xooii^) = Cuiii^) + 2; 

G2 : Xi = {x) = Cuii (x) + Cui2 {x) + 2, X2 = Xuj2 i^) = Cuj2 ix) + l. 

Since products of characters decompose into their sum, the recursive construction can proceed, 
but the polynomials will be different. Both approaches clearly can rightfully claim to be 
generalizaiions of classical Chebyshev polynomials. 

• Our approach to the derivation of multidimensional orthogonal polynomials hinges on the 
knowledge of appropriate recursion relations. The basic mathematical property underlying 
the existence of the recursion relation is the complete decomposability of products of the orbit 
functions. Numerous examples of the decompositions of products of orbit functions, involving 
also other Lie groups than SU(n), were shown elsewhere [15, 16]. An equivalent problem is 
the decomposition of products of Weyl group orbits [7] . 
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• The possibility to discretize the polynomials is a consequence of the known discretization of 
orbit functions. For orbit functions it is a simpler problem, in that it is carried out in the 
real Euclidean space M". In principle, it carries over to the polynomials. But variables of the 
polynomials happen to be on the maximal torus of the underlying Lie group. Only in the 
cases of (3.3), the variables are real (the imaginary unit multiplying the 5-functions can be 
normalized away). For An with n > 1, only some of the orbit functions of the congruence 
class are real valued. Practical aspects of discretization deserve to be further investigated. 

• For simplicity of formulation, we insisted throughout this paper that the underlying Lie 
group be simple. The extension to compact semisimple Lie groups and their Lie algebras is 
straightforward. Thus, orbit functions are products of orbit functions of simple constituents, 
and different types of orbit functions can be mixed. 

• Polynomials formed from £^-functions by the same substitution of variables should be equally 
interesting once n > 1. We know of no analogs of such polynomials in the standard theory 
of polynomials with more than one variable. Intuitively, they would be formed as sums 
C + S polynomials, Their domain of orthogonality is twice as large as that of Chebyshev 
polynomials. The iiJ-functions have been studied in [11, 17, 42]. 

• The generating functions for characters of irreducible finite dimensional representations of 
simple Lie groups were invented in [43]. They contain wealth of information about the 
characters and therefore also about the orbit functions and hence about the polynomials. In 
principle they arc not limited by the number of variables. Unfortunately, these functions 
grow rapidly in complexity for higher groups. The example of the generating function of G2 
characters is rather convincing [6]. However there is an aspect to the character generators 
that is not matched in other methods of constructing the polynomials. They provide rather 
specific information about existence of syzygies (identities) in the polynomial rings. 

APPENDIX 

A Orbit functions of Ai, their polynomial forms and Chebyshev polynomials 

A number of multivariate generalizations of classical Chebyshev polynomials are available in liter- 
ature [4, 27, 28, 44, 45], the aim of this section is to show in all details how Chebyshev polynomials 
appear as particular case of the multivariate polynomials proposed in this paper. First we recall 
that well-known classical Chebyshev polynomials can be obtained independently using only the 
properties of C- and S-orbit functions of the Lie group Ai, see [39] for detailes. The C-polynomials 
yielded by our approach are naturally normalized in a different way than the classical polynomials. 
It makes the correspondence between C- and S'-polynomials a direct special case of the general 
properties (2.4) and (2.5). 

We start a derivation of Ai polynomials from the beginning in a way that emphasizes the 
underlying Lie algebra Ai and, more importantly, in a way that directly generalizes to simple Lie 
algebras of any rank n and any type, resulting in polynomials of n variables and of a new type for 
each Lie algebra. 

The construction yields a different normalization of polynomials (form of Dickson polynomials) 
and their trigonometric variables than is common for classical Chebyshev polynomials. In ID, no 
new polynomials emerge than those equivalent to the classical Chebyshev polynomials of the first 
and second kind. Our insight into the structure of the problem for the general simple Lie group 
was gained in this construction. We are inclined to refer to the Chebyshev polynomials derived in 
this paper as the canonical ones. 
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The underlying Lie algebra Ai is often denoted su{2). In fact, this case is so simple that the 
presence of the Lie algebras never needed to be acknowledged^. 

The size of an orbit of Ai with the highest weight (m) in cj-basis is either 2, if m > 0, or 1, 
if m = 0. A weight belongs to the congruence class specified by the value #(m) = m mod 2. 
In particular, all exponents of monomials in a polynomial Cm or Sm have the parity of m. The 
dimension of the representation (m) is d^m) = m + 1. 

The orbit functions of Ai are of two types: 

C^(x) = e2™ + e-2~ = 2cos(27rmx), xeR, m G Z>°; (A.l) 
Sm{x) = e^''^"*^ - e-2~ = 2i sin(27rma;), x eR, m G Z>°. (A.2) 

When m = 0, the general definitions (2.2) and (2.3) give Co{x) = 1, and So{x) = for all x, 
but (A.l) gives Co{x) = 2. We keep this convention in this section. 

The simplest substitution of variables that would transform the orbit functions into a polyno- 
mial is Z"^ = e^'^*™^. Exponents in such a polynomial are the integers m and —m. Chebyshev 
polynomials can be represented in symmetric way: with symmetric positive and negative powers, 
and with all coefficients equal to 1 or —1. Instead, we introduce new variables X and S as follows: 

X := Ci(x) = e^™ + e-^'^'^ = 2cos(27r.x), 

S := Si{x) = e^"*^ - 6-2^^*^ = 2isin(27rx). ^ ^ 

Polynomials can now be constructed recursively in the degrees of X and S by calculating the 
decompositions of products of appropriate orbit functions. 'Generic' recursion relations are those 
where one of the first degree polynomials, X or S, multiplies the generic polynomial Cm or Sm, i-e. 
m > 1. Omitting the dependence on x from the symbols, we have the generic recursion relations 

^Cm = Cm+1 + Cm-1, SCm = Sm+l — Sm-1, m > 1 

XSm = Sm+l + Sm-1, SSm = Cm+1 — Cm—1, 

When solving recursion relations for C-polynomials, we need to start from the lowest ones: 

X^ = C2 + Co = C2 + 2 =^ C2 = X^-2, 

XC2 = C3 + Ci = C3 + X =^ C3 = XC2 — X = X"^ — 3X, 

XCm = Cm+1 + Cm-1 ^ Cm+1 = XCm ~ Cm-1, m ^ 3. 

Several lowest results are in Table 3. Hence we conclude that Cm = 2Tm, for m = 0, 1, . . . . 

There are also recursion relations for C-polynomials resulting from products of two S-polyno- 
mials: 

S^ = C2-2 ^ C2 = S^ + 2, 

SSm = Cm+1 — Cm—1 =^ Cm+1 = SSm + Cm—1, 2. 

Note that each Cm can be written in two ways, as a polynomial of degree m in X, and as 
a polynomial of the same degree involving S and X. Equating the two expressions for Cm, we 
obtain a trigonometric identity for each m. For example, from the two ways of writing C2, we find 

- 5-2 = 4 sin2(27rx) + cos^{2'kx) = 1. 

Our S is defined to be purely imaginary, hence the negative sign at S"^. 



^The Lie algebra su{2) is the backbone of the angular momentum theory in quantum physics. The normalization 
of its generator implied here is common in mathematics. In particular, m equals to twice the angular momentum in 
physics. 
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X'^ - 6X^ + lOX^ - AX 


Table 3: The irreducible C- and S'-polynomials of of degrees up to 8. 



Analogous relations yield polynomial expressions for Sj, 



XS = S2 

XS2 = 83 + 8 

X8m = Sm+1 + Sm-1 

SC2 = 83 — 8 

SCm = Sm+\ + Sm—1 



82 = XS 

8"^ = X82 — S = X^ S — S 

8m+i = 8Sm — 8fn—i, m ^ 2; 

Ss = SC2 + S = S^ + 38 

8m+l = 8Sm — Cm-1, Ul ^ 2. 



A fundamental relation between 5*- and C-polynomials (appropriately normalized) is the special 
case of (2.4). It eliminates the need to find S'-polynomials recursively, provided the C-polynomials 
have been found. The character Xmix) of an irreducible representation of Ai of dimension m + 1 
is known explicitly for all m ^ 0. There are two ways to write the character: as the ratio of 
S-functions, and as the sum of C-functions. Explicitly, that is 



Xni{x) 



Cm{x) + Cm-2{x) + ■ ■ ■ + 



C2{x) + Co if m even. 



8{x) ■ ■ ■ [C3{x) + X{x) if m odd. 

Note that (A. 4) is the Chebyshev polynomial of the second kind Um{x). 



(A.4) 



Remark 4. Note that in the character formula we used Cq = 1, while for C-polynomials we used 
Co = 2. It is just a question of normalization of orbit function Cq. Here we used it in order to 
obtain classical form of the Chebyshev polynomials. More generally for any simple G, it is sometime 
convenient to scale up orbit functions of non-generic point, say A, by the factor equal to the order 
of the stabilizer of A in the Weyl group W. 



Remark 5. The main argument in favor of our normalization of Chebyshev polynomials is that 
polynomials Cm from Table 3 are Dickson polynomials (it is well known that them are equivalent 
to Chebyshev polynomials over the complex numbers). It is easy to prove (see e.g. [39]) that Weyl 
group of An is equivalent to Sn+i, therefore it is natural to consider multivariate C-polynomials 
of An as n-dimensional generalizations of Dickson polynomials (as permutation polynomials) . Also 
our form of Dickson-Chebyshev polynomials makes them the lowest special case of (2.4) and (2.5) 
without additional adjustments and it appears more 'natural' because, for example, the equality 
C| = C4 -|- 2 would not hold for T2 and T4. 
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B Recursion relations for A2 orbit functions and polynomials 

In the previous section, variables X and S played an almost symmetrical role. This is not the case 
when the rank of the Lie algebra L„ exceeds 1. The number of variables of type X is equal to n, 
the number of exponential functions comprising such variables is a (small) divisor of the order \W\ 
of the Weyl group. In contrast, the variable S is unique for all It is a sum of the maximal 
number of exponential functions, namely \W\. 

The variables of the A2 polynomials are the C-functions of the lowest dominant weights cji = 
(1, 0), 0J2 = (0, 1), and the unique lowest S'-function whose dominant weight is (1, 1). The variables 
are denoted as follows, 

Xi := C(^i^Q){xi,X2), X2 := C^o,i){xi,X2), S := S'(i,i)(xi, ^2). (B.l) 

We omit writing (xi,X2) at the symbols of orbit functions for simplicity of notations. 

In addition to the obvious polynomials Xi, X2, Xf, X1X2, and X2, we recursively find the rest 
of the ^2-polynomials. The degree of the polynomial C(a 5) equals a + b. The degree of S(^a.,b) is also 
a + b provided ab ^ 0, otherwise the ^-polynomials are zero. 

Due to the A2 outer automorphism, polynomials C(a,6) and C(5 are related by the interchange 
of variables Xi -H- X2 (i.e. C(„,6)(Xi, X2) = C(6^„)(X2, Xi)). 

In general, each term in an irreducible polynomial, cquivalcntly each weight of an orbit, must 
belong to the same congruence class specified by the congruence number #. For A2-weight (a, 6), 
we have 

#(0, b) = {a + 2b) mod 3. (B.2) 

Hence, irreducible orbit functions have a well defined value of ^. For yl2-orbit functions, we have 
#(C((j 5)) = #(5'(a,6)) = (a + 26) mod 3. Consequently, there are three classes of polynomials 
corresponding to # = 0,1,2. During multiplication, the congruence numbers add up mod 3. 
A product of irreducible orbits decomposes into the sum of orbits belonging to the same congruence 
class. 

The sizes of the irreducible orbits of ^(^2) are found in Table 2. The dimension d^a,b) of the 
representation of A2 with the highest weight (a, b) is given by 

d^a,b) = hi^ + '^){b + l){a + b + 2). (B.3) 



B.l Recursion relations for C-function polynomials of A2 

There are two 4-term generic recursion relations for C-functions. They arc obtained as the de- 
composition of the products of X and Y, each being a sum of three exponential functions, with 
a generic C-function which is the sum of |W^(^2)| = 6 exponential terms. We call C(„ ;,) generic, 
provided it is a sum of 6 distinct exponential terms. In order that a recursion relation be generic, 
the product has to decompose into three distinct generic (i.e. 6-term) C-functions 

XlC(a,b) = C(a+l,b) + C'(a-1,6+1) + Ci^a,b-1), a, 6 > 2; 

-'^2C(a,6) = C'(a,6+1) + C'(a+1,6-1) + C'(o-l,6)) a,b^2. 

Before generic recursion relations can be used, the special recursion relations for particular values 
a, 6 G {0, 1} need to be solved recursively starting from the lowest ones: 

Xf = C(2,o) + 2X2, X2 = C(o,2) + 2X1, X1X2 = C(i^i) -I- 3; 

-^^1^(1,1) = C'(2,i) + 2C(o,2) + 2X1, X2C(i^i) = C(i_2) + 2C(2,o) + 2-^2; 

for a ^ 2 : 
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-^^1^(0,1) = ^"(0+1,1) + C'(a-1,2) + 2C(„^o)) -^^2(^(0,1) = C'(a,2) + 2C(„+i^o) + (^(a-l,!); 

-'^lC'(a,0) = ^"(0+1,0) + C'(a-l,l)) -^^2(^(0,0) = C'(a,l) + C(a-1,0); 

for 6 ^ 2 : 

-^^1^(1,6) = C(2,b) + 2C(o,6+i) + C(i^ft_i), X2C(i^ft) = C(i^ft+i) + C(2,b_i) + 2C(o,6); 

-^lC'(0,fe) = ^"(1,6) + (^(0,6-1), X2C(o,6) = (^(0,6+1) + C'(l,6-1)- 

Using the symmetry of orbit functions with respect to the permutation of the components of 
dominant weights, we obtain analogous polynomials C(^afl) ^^^d 0(^^,1) for a G N. Then the 4- 
term special recursion relations are solved yielding C(^2,b) ^■^d (^(0,2) a,b eN. After that, the 

generic recursion relations should be used. 

B.2 Recursion relations for S'-function polynomials of A2 

Polynomials >S'(a^o) = '^'(0,6) = vanish for all a and b. There are two 4-term generic recursion 
relations for S'-functions, 

XlS(^a,b) = S(a+l,b) + S(^a-l,b+l) + S(a,b-1)^ a, 6 ^ 2; 

^2S(^a,b) = S{a,b+1) + 'S'(o+l,6-l) + S(^a-l,b)j a, 6 ^ 2. 

Then there are the special 3- and 4-term recursion relations for particular values of a = 1 and/or 
6 = 1: 

XiS = 3(^2,1), = 3(^1^2)', 

-^l'S'(l,2) = '5(2,2) + S, X23{i^2) = '5(1,3) + '5(2,1); 

-'^l '5(2,1) = '5(3,1) + '5(1,2), -^^2.5(2,1) = 5(2,2) + 

-'^l -5(3,1) = -5(4,1) + -5(2,2), -^^2-5(3,1) = 3^2,1) + -5(3,2)- 

Further on, generic recursion relations can be used. 

There are additional recursion relations for C-polynomials, even if not necessarily used for finding 
the polynomials: 

S3 = C{2,2) - 2C(o,3) - 2C(3^o) + 2C(i^i) - 6, 

'5'5(2,l) = C(3^2) - (^(1,3) - 2C(4^o) + 2C(o,2) - 2C(i^o) + C'(2,l), 
'5'5(i,2) = (^(2,3) - C(3^i) - 2C(o,4) + 2C(2,0) - 2C(o,l) + C{i^2)-, 
S3(2,2) = C(3,3) - (^(4,1) - ^1,4) + 2C(o,3) + 2C(3^o) " '^(l,!); 

for a, 6 ^ 3 : 

S3(^a,b) = C'(a+l,b+l) — C(^a+2,b-l) — C'(a-l,6+2) + ^"(0-2,6+1) + C'(a+i,b_2) — C(a_l,6_l). 

These are interesting for a different reason. Indeed, a particular 3(a,b) can be written as a linear 
combination of C-functions divided by 3. From (2.4), we find that 3(^a,b) is written as a polynomial 
of C-functions multiplied by 3. 

B.3 The character of A2 

The character X(a,b) is given either as a fraction of S'-functions (Weyl character formula) or as 
a linear combination of the C-function. In the A2 case the general formula (2.4) is specialized 

X{a,b) {x, y) = '^(°+^'''+^)^^'^^ = C(„ b) {x, y) + V mxCx{x, y). 
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C-polynomials 


# = 


C(i,i) 


- 3 


C(3,0) 


- 3X1X2 + 3 


C(2,2) 


X^Xl - 2X1 - 2^1 + 4X1X2 - 3 


# = 1 


'^(1,0) 


Xi 


C(0,2) 


X2 — 2X\ 


^"'(2.1) 


Xf X2 - 2X| - Xi 


C(l,3) 


XiX| - 3x^X2 - x| + 5X1 


<^(4,0) 


Xf - 4X^X2 + 2X| + 4X1 



Table 4: The irreducible C- and S'-polynomials of 
or S(^a,b) we obtain C(b,a) or S'(b,a) respectively by 



S'-polynomials 


# = 


•^(2,2) 


X1X2 - 1 


5'(1,4) 


X| - 2X1X2 + 1 


'S'(3,3) 


y2 v2 vZ vZ 
^1^2 ~ ^1 ~ ^2 


# = 1 


•5(2,1) 


Xi 




x| - Xi 


-5(3,2) 


XfX2 - xl - Xi 


•5(2,4) 


X1X2 — 2X^X2 — x'2 "l~ 2X2 


-5(5,1) 


x^ — 3x^X2 + X| + Xl + X2 



A2 of degree up to 4. Prom any polynomial C((j 
interchanging Xi and X2. 



The summation extends over the dominant weights that have positive multiplicities m\ in the case 
of X(a,6)- The coefficients (dominant weight multiplicities) are tabulated in [2] for the 50 first X(a,6) 
in each congruence class of A2. The first few characters for the congruence class # = are: 



X(o,o) 


= ^"(0,0) 


= 1, 


X(i,i) 


= C(l,l) 


+ 2C(o,o) , 


X(3,0) 


= C'(3,0) 


+ C'(ia) + C'(o,o)) 


X(0,3) 


= ^"(0,3) 


+ C'(i,i) + C'(o,o)) 


X(2,2) 


= ^"(2,2) 


+ (^(0,3) + C'(3,o) + 2C(i_i) + 3C(o,o), 


X(l,4) 


= C(i,4) 


+ C'(2,2)) + 2C(o,3) + (^(3^0) + 2C(i^i) + 2(7(0,0), 


X(3,3) 


= ^"(3,3) 


+ C(4,i) + C(i,4) + 2C(2,2)) + 2C(o,3) + 2C(3,o) + 3C(i_i) + 4C(o,o) 


X(6,0) 


= <^(6,0) 


+ C'(4,l) + C'(2,2)) + C'(0,3) + C'(3,0) + (^(l.l) + C'(0,0)- 



The equalities must satisfy two relatively simple conditions: (i) The dominant weights on both 
sides must have the same congruence number (B.2), and (ii) the number of exponential terms 
in a character X(a,6) is known to be the dimension (B.3) of the irreducible representation (a, 6). 
Therefore, the sizes of the orbit functions on the right side have to add up to the dimension. For 

# = 1: 

X(i,o) = ^"(1,0), 

X(o,2) = ^"(0,2) + C'ci.o), 

X(2,l) = C{2X) + C'(0,2) + 2C(i,o), 

X(l,3) = ^"(1,3) + C(2,l) + 2C(o,2) + 2C(i,o), 

X(4,0) = ^"(4,0) + C'(2,l) + ^"(0,2) + C(1,0), 

X(0,5) = ^"(0,5) + C'(l,3) + ^2,1) + C'(0,2) + C'(1,0)- 

For # = 2, it suffices to interchange the component of all dominant weights in the equalities for 

# = 1. Thus no independent calculation is needed. 

Example 2. Comparison of polynomials obtained by our exponential substitution method in the 
case of A2, with the results of [19] (III), reveals coincidence of the polynomials in both cases, as 
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demonstrated in this example. Moreover, the polynomials considered here are the partial case of 
the Mecdonald polynomials. 

Suppose X = {xi,X2) is given in the a-basis, then C-functions A2 of assume the form: 



C(ofi){.x) = 1, 

C{o,m){x) = C(^,o)(x) = e-2~i + e2~ie-2~2 + e^™^, 

( \ 2nimixi 2nim2X2 i —2nimixi 2iTi{mi+m2)x2 i 

'-'{mi,m2)\-^} ~^ ^ "T ^ ^ "T 

^27ri{mi+m2)xi ^—2iTim2X2 _j_ ^—2nim2Xi ^—2TrimiX2 _j_ 
g— 27ri(mi+m2)a;ig27rimia;2 _j_ ^2irim2Xi ^—2TTi{mi+m2)x2 



^27rimixi ^2Trim2X2 ^—2Trimixi ^2Tri(mi+m2)x2 

g27ri(mi+m2)j;ig-27r-tm2X'2 _ ^-2TTim2Xi ^-2T7im-i_X2 _j_ 
g— 27ri(mi+m2)a;ig27rimia;2 _j_ g27rim2a;ig— 2iri(mi+m2)a;2 



The polynomials e''" and e~ given in (2.6) of [19, III] coincide with those in (B.4) whenever the 
correspondence a = 27ra;i, r = 27r.T2 is set up. Thus both our orbit function polynomials of A2 and 
of [19, III] are orthogonal on the interior of Steiner's hypocycloid. 

It is noteworthy that the regular tessellation of the plane by equilateral triangles considered 
in [19] is the standard tiling of the weight lattice of A2. The fundamental region R of [19] coincides 
with the fundamental region F(^A2) in our notations. The corresponding isometry group is the 
affine Weyl group of A2. 

Furthermore, continuing the comparison with [19], we should emphasize that it is known that 
orbit functions are eigenfunctions not only of the Laplace operator but also of the differential 
operators built from the elementary symmetric polynomials, see [15, 16]. 

C Recursion relations for C2 orbit functions 

There are two congruence classes of C2 orbit functions/polynomials. For C2 weight (a, b) (dominant 
or not), we have 

i^{a,b) = a mod 2 (C.l) 
The dimension c?(a,6) of an irreducible representation of C2 with the highest weight (a, b) is given by 
d^a,b) = U^ + 'i-){b+l){2a + b + 3){a + b + 2). (C.2) 



In multiplying the polynomials, congruence numbers add up mod 2. Character in the case of 
C2 is given by (2.4), where the C- and S-functions are those of C2, as are the coefficients m\ (also 
tabulated in [2]). 

We denote the variables of the C2-polynomials by 

Xi := C(i^Q){x,y), X2 := C(o,i)(a;, y), and S := S'(i_i)(a;, y), 

often omitting (x, y) from the symbols. The variable S cannot be built out of Xi and X2. Although 
the variables are denoted by the same symbols as in the case of A2 (and also G2 below), they are very 
different. Thus Xi and X2 contain 4 exponential terms and S contains 8 terms. The congruence 
number # of Xi and 5 is 1, while that of X2 is 0. 
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C-polynomials 


# = 


^(0,1) 




C(2,0) 


X? - 2X9 - 4 


^(2,1) 


y\ -1 /\ ^ y\ i-\ /A 


C(4,0) 


4 - 4X2 + + 8X2 - 4X^X2 + 2X| 


C(0,2) 


4 - 2X2 ^ ^ ^2 


C(0,3) 


9X2 - 3XfX2 + 6X| + X| 


C(2,2) 


-8 + 10X2 - 2Xf - 20X2 + 8X2X2 - 12X2 ^ xlXl - 2X| 


C(0,4) 


4 - 8X2 _^ 2X1^ + I6X2 - 8X2X2 + 20X1 - 4X2x1 + 8X| + X| 


# = 1 


C(i,o) 


^1 




X1X2 — 2X1 


C(3,0) 


Xj' — 3X]^X2 — 3X]^ 


<^(3,1) 


2X1 - 4X1X2 + Xi'X2 - 3XiX| 


C(l,2) 


6X1 - 2Xf + 3X1X2 + XiX| 


C(l,3) 


-6X1 + 2Xf + 6X1X2 - 3Xf X2 + 5X1 X2 + XiX| 



Table 5: The irreducible C-polynomials of C2 of degree up to 4. 



C.l Recursion relations for C-functions of C2 

The two generic recursion relations for C-functions of C2 are 

XlC(^a,b) = C{a+l,h) + C{a-l,b+l) + C{a+l,b~l) + C{a-l,b)i O, 6^2; 

X2C(a,b) = C{a,b+l) + C'(a+2,6-l) + C'(a-2,6+l) + (^(a.ft-l)) a ^ 3, 5^2. 

The special recursion relations for C-functions involving low values of a and b have to be solved 
first starting from the lowest ones: 



-^^1(^(1,6; 

X1X2 = 

-^2C(aj 
-'^2C(a^0 

-^2C(o,6; 
X2C(2,1 
X2C(i^l 

-'^2C(2,0 



— C(^a+l,l) + ^0-1,2) + 2C((j+i^o) + C(^a-l,l)i 
= C{^a+lfi) + C{a-l,l) + C(a-1,0)) 



C(2,b) + 2C| 
C(i,b) + C 



(0,6+1) 



+ ^"(2,6-1) + 2C, 



(0,6) > 



(1,6-1), 

= C(2,i) + 2C(o,2) + 2C(2,o) + +2X2; 
C(i,i) + 2X1; X2 = C(2,o) + 2X2 + 4; 

= C{^a,2) + 2C(a+2,0) + ^"(0-2,2) + 2C(a^o)> 



C, 



(a,l) 



(a-2,1)' 



C(2,6+l) + C(4^fe 
C(l,6+1) + ^"(3,6 
C(0,6+l) + ^^(2,6-1) + C'l 



1) + 2C(o,6+i) + C(2,fe_i), 
1) + (^(1,6-1) + C'(l,6), 

(0,6-1), 

C(2,2) + 2C(4^o) + C'(0,4) + 2C(o,2) + 2C(2,o); 

C(i,2) + 2C(3^o) + C'(i,i) + 2X1; 



C(2,l) + 2X2; 



xl 



C(o,2) + 2C(2,o) + 4- 



a ^2; 

a ^ 2; 

b^2: 



o ^ 3; 
o ^ 3; 
b^2; 
6^2; 
6^2; 



The 3- and 4-term recursion relations are solved independently, giving us C(o,(,), (^(0,0)! ^(1,6)5 
and C((j 1) for all a and b, e.g. see Table 5. 
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C.2 Recursion relations for S^-functions of C2 

The generic relations for S-functions are readily obtained from those of C-functions by replacing C 
by S, and by making appropriate sign changes, 



ia+l,b) 



S{a-l,b+l) + 'S'(a+l,fe-l) — 



^2S(^a,b) — '5(0,6+1) - '5(0+2,6-1) + 'S'(a-2,6+l) " 'S'(a,6-1)) 



a, 6^2; 
a > 3, b^2. 



The special recursion relations for S'-functions involving low values of a and b have to be solved 
first, starting from the lowest ones: 



'5(0+1,1) + •5(o-l,2) + •5(o-l,l)) 



-'^i'5(o,i) 

-'^2'S'(a,l) = -5(0,2) + -5(0-2,2)) 

^2-5(2,6) 
-^^1-5(1,6) 
-^^2-5(1,6) 



-5(2,6+1) + -5(4,6-1) + -5(2,6-1); 

-5(2,6) + -5(2,6-1)) 

-5(1,6+1) + -5(3,6-1) + -5(1,6-1) - •5(1,6)) 
-^^2-5(2,1) = -5(2,2); 

XiS = -5(2,1); 

X2S 



a > 2 
a ^ 3 
6^2 
5^2 
b>2 



5(1,2) - S. 



SS^ 



(0,6) 



(a+1,6+1) 



(a-1,6+2) 



C(a+3,6-l) — C(^a-3,b+l) 



+ C'(o-l,6-l) - C'(a+l,6-2) + ^^(0+3,6-2) + C'(o-3,6+2)) 



a > 4, 6^3; 



SS = C(2,2) - 2C(o,3) 



2C, 



(4,0) 



2^^(2,0) ~ 2X2 + 2X2 + 2C(-2,i) + 8. 



All C-functions of C2 are real-valued. Here are a few examples of C2-characters: 



# = 0: 

X(o,o) = 

X(o,i) = 

X(2,0) = 

X(0,2) = 

X(2,l) = 

X(0,3) = 

X(4,0) = 

X(2,2) = 

# = 1: 

X(i,o) = 

X(i,i) = 

X(3,0) = 

X(l,2) = 

X(3,l) = 

X(l,3) = 



C(o,o) = 1; 

1 + ^"(0,1) = 1 + X2; 

2 + X2 + C(2,o); 

2 + X2 + C(2,0) + (^(0,2); 

3 -I- 3X2 -I- 2C(2,o) + ^"(0,2) + ^"(2,1); 

2 -I- 2X2 -I- C(2,o) + (^(0,2) + ^"(2,1) + (^(0,3); 

3 -I- 2X2 -I- 2C(2,o) + ^"(0,2) + ^"(2,1) + (^(4,0); 

5 -I- 4X2 -I- 4C(2,0) + 3C(o,2) + 2C(2,i) -I- C(o,3) + C(4^o) + C'(2,2)- 



(1,0) 



^1; 



2X1 + C(i,i); 

2X1 + C(i,i) + C(3,o); 

3X1 + 2C(i,i) -I- C(3,o) + (^(1,2); 

4X1 + 3C(i,i) + 2C(3^o) + Ql,2) + ^3,1); 

4X1 + 3C(i,i) -I- 2C(3,o) + 2C(i,2) + C(3,i) C(i,3). 



Using these characters and Table 5, we can calculate all irreducible S'-polynomials of degree up 
to four with respect to the variables Xi and X2 using (2.5). Note that X(o,4) yields the polynomial 
of order five. 
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S'-polynomials 


# = 


^(1,2) 


I+X2 


'^(3,1) 


-2 + X? - Xo 


-^(1,3) 


2-Xf + 3X2 + Xl 


■^(3,2) 


— 1 — + X-^ X2 — X2 


-5(1,4) 


2-Xl + 7X2 - 2x1X2 + ^ xl 


5(5,1) 


3 - 4X^2 + Xf + 4X2 - 3X2^2 + X22 


•5(3,3) 


-3 + 4X2 _ xf - 7X2 + 3x2X2 - 5X2 ^ xlXl - X| 


# = 1 


-5(2,1) 


Xi 


-5(2,2) 


X1X2 


-5(4,1) 


— 3X2 ~t~ x^ — 2x^^X2 


-5(2,3) 


2X1 — X^ -\- 2X1X2 ~\- X1X2 


-5(4,2) 


-4X1X2 + XfX2 - 2X1X22 


-5(2,4) 


5X1X2 - 2Xf X2 + 4X1 X| + XiXf 



Table 6: The irreducible S'-polynomials of C2 of degree up to 4. 



D Recursion relations for G2 orbit functions 

All G2 weights fall into the same congruence class # = 0. Thus there are no congruence classes to 
distinguish in G2. The dimension of the irreducible representation (a, 6) of G2 is given by 

d{a,b) = Tio (« + + + ^ + 2)(2a + & + 3)(3a + b + 4)(3a + 2b + 5). (D.l) 

The variables are the orbit functions of the two fundamental weights coi = (1,0) and UJ2 = (0, 1). 
Also S{x,y) := S'(i i)(x,?/) serves as an independent variable 

^1 ■= C(i^o){x,y), X2 := C(o,i)(a;,y), S = S'(i,i)(x, y). 
The variables Xi and X2 are the sums of 6 exponential terms, while S has 12 terms, also 

D.l Recursion relations for C-functions of G2 

There are two generic recursion relations for C-polynomials of G2, each containing one product 
term and six C-polynomials: 

for a^3, 6^4 : 

^lC(a,b) = C'(a+1,6) + C'(a-l,6+3) + ^"(0+2,6-3) + C'(a-2,6+3) + C'(a+l,6-3) + ^"(0-1,6); 

for a^2, 6^3 : 

X2C(^a,b) = C'(a,b+1) + ^"(0+1,6-1) + ^"(0-1,6+2) + C(^a+l,b-2) + C'(a-1,6+1) + C'(a,6-1)- 

Specializing the first of the generic relations to either a G {0, 1, 2} or 6 G {0, 1, 2, 3}, we have 

-^^1^(2,6) = C'(3,6) + C'(l,fe+3) + ^^(4,6-3) + 2C(o,6+3) + C(^3^h-3) + ^'(1,6)! 
-^^1^(1,6) = C'(2,b) + 2C(o,b+3) + C(^3^b-3) + C(2,b-3) + (^(1,6) + 2C(o,b); 
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-^^1^(0,6) = ^7(1,6) + C'(2,&-3) + C'(l,6-3); 

-^^1^(0,3) = C'(a+i,3) + C{a-l,e) + 2C(a+2,0) + ^"(0-2,6) + 2C(a+l,0) + C'(a-1,3); 

-^^1^(0,2) = ^"(0+1,2) + + C'(a+l,l) + ^"(0-2,5) + C'(a,l) + (^(0-1,2); 

-^^1^(0,1) = C'(a+l,l) + C'(a-1,4) + ^"(0-2,4) + C'(a-l,l) + C'(a,2) + (^(0-1,2); 

-'^lC'(a,0) = ^"(0+1,0) + C'(a-1,3) + ^"(0-2,3) + C'(a_i,o); 

-^^1^(2,3) = C'(3,3) + C'(l,6) + 2C(4^o) + 2C(o,6) + 2C(3_o) + (^(1,3); 

-^^^1(^(2,2) = (^(3,2) + (^(1,5) + (^(0,5) + C'(l,2) + (^(2,1) + C^s^i)] 

XlC(^2,l) = C'(3,l) + C'(l,4) + 2C(o,4) + C'(l,l) + (^(2,2) + C'(l,2); 

-^^1^(2,0) = C'(3,o) + C(i^3) + 2C(o,3) + Xi; 

-^^1^(1,3) = ^"(2,3) + 2C(o,6) + 2C(3^o) + 2C(2,0) + ^"(1,3) + 2C(o,3); 
-^^1^(1,2) = ^"(2,2) + 2C(o,5) + C'(l,2) + 2C(o,2) + C'(2,l) + C'(l,l); 

-^^1^(1,1) = C(2,i) + 2C(o,4) + ^"(1,2) + 2C(o,2) + + 2^2; 

^lC{0,3) = ^"(1,3) + 2C(2,o) + 2^1; 
-^^1^(0,2) = C'(l,2) + ^"(1,1) + 2-^2; 

XiXi = C(2,o) + 2C(o,3) + 2X1 + 6; 

X1X2 = (^(l,!) + 2C(o,2) + 2X2. 

Specializing the second of the generic relations to either a G {0, 1} or 6 G {0, 1, 2}, we have 

X2C(^l^b) = C'(l,6+l) + ^"(2,6-1) + 2C(o,6+2) + C'(2,6-2) + 2C(o,6+i) + C(i_6_i); 
-'^2C'(0,6) = ^"(0,6+1) + C{l,b-1) + C'(l,6-2) + C{0,b-iyi 

X2C(^a,l) = C(a,2) + 2C(a+i^o) + C'(a-1,3) + C'(a-1,2) + 2C(a,0) + C'(a,l); 
X2C(^a,0) = C(a,l) + ^"(0-1,2) + C(^a~l,l)', 

X2C(^l^2) = C(i^3) + C(2,l) + 2C(o,4) + 2C(o,3) + C(i^i) + 2C(2,o); 

X2C(^i^i) = C(i_2) + 2C(2,o) + 2C(o,3) + 2C(o,2) + C'(i,i) + '^Xi; 

X2C(^0,2) = ^"(0,3) + C'(l,l) + 2-^^! + X2; 

X2X2 = C{o,2) + ^Xi + 2X2 + 6. 
Remark 6. It can be seen from Table 7 that order of C(a (,)-polynomial sometimes exceeds a + 6. 
D.2 Recursion relations for S'-functions of G2 

Generic recursion relations for S'-polynomials differ very little from those for C-polynomials. 

for a^3, 6^4 : 

XlS(^a,b) = S(a+l,b) + S(a-l,b+3) + 'S'(a+2,6-3) + •S'(a-2,6+3) + •S'(a+l,b-3) + 5'(a-l,()); 
for a^2, 6^3 : 

X2S(^a,b) = 'S'(a,6+1) + 'S'(a+l,b-l) + 'S'(a-l,6+2) + 'S'(a+l,fe-2) + 'S'(a-l,b+l) + 'S'(a,b-1)- 

The 5-polynomials need not be calculated independently. They can be read off the tables [2] as 
the characters of G2 representations. The dimension d(^a,b) of irreducible representation of G2 
with the highest weight A = {a,b) is given by (C.2). 

Here are all G2-characters X{a,b) with a + 6 ^ 3: 

X(i,o) = 1 + Qi.o) = 1 + -'^^i; 



22 



C-polynomials 








X2 




-6 - 2X1 - 2X2 + Xl 




12 + 4X1 + 2X2 + X1X2 - 2X1 




-12 - 10X1 - 2X^2 _ ^X2 - 3X1X2 + 2X| + XiX| 




-12 - 12X1 - 2X2 _ 3^ _ 3X1X2 + 2X2 _^ XiX"^ 


C(2,0) 


18 + 22X1 + 5X^ + 6X2 + 6X1X2 - 4X| - 2X1 X| 


C(l,3) 


-36 - 58Xi - 22X2 _ ~ I2X2 - I5X1X2 - 3x2X2 + 8X3 + 6X1X3 + XiX| 


^(0,4) 


6 + 8X1 + 2Xi2 - 8X9 - 10X1 X9 - 2Xi2X9 - 4X^ - 4X1 X^ + 2X^ + Xi X^ 


C(2,l) 


6X1 + 2X^2 _^ 20X2 + 24X1X2 + 5X^X2 + 6X| + 5X1 X| - 4X| - 2X1 X| 


(^(1,4) 


-4X^2^1 + 4X| + 4XiX| + Xi2x| 


'-'(3,0) 


fin -S^OiQY, -1- 48 Y'^ 7 )r3 i i c y„ i 97 v _i_ q \^2 1 9 )f 2 1 9 )r, )r2 t^2 v-2 


C(2,2) 


-108 - 156Xi - 46X^2 - 2Xf - 64X2 - 6OX1X2 + 9X^2X2 + 6X^X2 + 32X| + 36XiX| 
+11X^2 X2 + 12X| + 6XiX| - 2X2 X| _ 4^4 _ 2XiX| 


^(3,1) 


ino 1 ^ ^cw i A A i o\^3 i a c\a v i a v v i o A v i o\^3'v^ on v2 
lUo + ioUAi + 44A]^ + ZX^ + iU4A2 + ioOAi A2 + o4A-^ A2 + ZAj^A2 — 2UA2 

-14XiX| - 2X2X1 - 20X| - 16XiX| - X2x| + 4X| + 2XiX| 


C(4,0) 


198 + 400X1 + 282X^2 + 84Xf + 9Xf + 240X2 + 36OX1X2 + 152X^2X2 + 20X1^X2 
+26X| + 28XiX2 - 8X2X1 - 4X^X1 - 52X| - 46XiX| - 8X^X1 
-8X| - 8X1 X| + 4X| + 2X1 X| 


Table 7: The irreducible C(a ;,)-polynomials of G2 with a + 6 ^ 4. 



X(o,i) = 2 + C(i^o) + C'cci) = 2 + Xi + X2; 

X(2,o) = 3 + 2X1 + X2 + C(2,o); 

X(i,i) = 4 + 4X1 + 2X2 + 2^(2,0) + (^(1,1); 

X(3,o) = 5 + 4X1 + 3X2 + 2C(2^o) + C(3,o); 

X(0,2) = 5 + 3X1 + 3X2 + 2C(2^o) + C'(i,i) + C'(3,0) + ^"(0,2); 

X(2,i) = 9 + 8X1 + 6X2 + 5C(2,o) + 3C(i,i) + 2C(3,o) + C(o,2) + C(2,i); 
X(i,2) = 10 + 10X1 + 7X2 + 7C(2,o) + 5C(i,i) + 3C(3,o) + 3C(3,o) + 2C(o,2) + 

2C(2,l) + C(4^o) + ^1,2); 

X(o,3) = 9 + 7X1 + 7X2 + 5C(2,o) + 4C(i,i) + 4C(3,o) + 3C(o,2) + 2C(i,i) + 

<^(4,0) + <^(1,2) + C'(3,l) + <^(0,3)- 



E Recursion relations for Lie algebras of rank 3 

E.l Recursion relations for C-functions of A3 

There are 4 congruence classes of A^, defined by 

#(a,6,c) = a + 26 + 3c mod 4. (E.l) 
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S'-polynomials 




l + Xi 




2 + X1+X2 


•5(3,1) 


21 + 24Xi + 5Xf + 7X2 + 6X1X2 - 4X| - 2X1 X| 


•5(2,2) 


52 -1- 52X1 -1- ^QX? -i- ifiXo -1- Xo — ^ox^ — 4X^ x^ 


•5(4,1) 


113 + 151X1 + 58X^2 + 7Xf + 35X2 + 40X1X2 + 9X1^X2 - 22X| - 16XiX| - 3X^2x1 


•5(1,3) 


107 + 148Xi + 58Xf + 7Xf + 33X2 + 4OX1X2 + 9X^X2 - 21X| - 16XiX| - 3Xf X| 


•5(3,2) 


249 + 332X1 + 123X? + 14Xf + 96X2 + IIIX1X2 + 23Xf X2 - 43X| - 29XiX| 
-6X^2 X| - 4X| - 2X1 X| 


•5(2,3) 


550 + 879Xi + 463X2 + 105Xf + 9Xf + 385X2 + 533XiX2 + 189X^X2 + 20X^X2 
-llXiX| - 32X| - 17X2 X| - 4X3 X| _ gQj^s _ soXiXf - 8X^X1 - 8X| 
-8X1 X| + 4X| + 2X1 X| 


•5(4,1) 


651 + 1063Xi + 543X^2 + 114Xf + 9Xf + 488X2 + 679XiX2 + 232Xf X2 
+22X^X2 - 32X1 X| - 51X| - 22X^X| - 4Xf X| - 80X| - 66X1 X| - 9Xf X| 
-4X| - 6X1 X| + 4X| + 2X1 X| 



Table 8: The irreducible -polynomials of G2 with a + 6 ^ 5. 



The dimension d(^a,,b,c) of the irreducible representation of with the highest weight A = (a, b, c) 
is given by 

d(a,b,c) = ^{a + l){b + l){c + l){a + b + 2){b + c + 2){a + b + c + 3). (E.2) 

The variables of the polynomials are chosen to be 

^1 ■= C(^ififi){xi,X2,X3), X2 := C(o,i,o)(a;i,a;2,X3), X3 := C(o,o,i)(a:i,a;2,a:3)- 

The orbit functions C(i 0,0) and (^(0,0,1) contain 4 exponential terms, and (^(0,1,0) has 6 terms. 
All terms have the form e'^'^^^^'^\ where fi runs over the points/weights of the corresponding orbit. 

As for C-functions, generic recursion relations are the decompositions of the following products, 
where we assume a, 6, c ^ 2: 

^lC(a,b,c) = C'(a+i,fe,c) + C{a-l,b+l,c) + C'(a,6-l,c+l) + C'(a,6,c-1); 
X2C^a,b,c) = C(^a,b+l,c) + C'(a+l,fe-l,c+l) + C'(a-l,6,c-|-l) + C'(o+i,;,,c-l) + 

C(a-l,6+l,c-l) + C'(a,6-l,c); 
X3C^a,b,c) ~ Qa,6,c+1) + C'(a,6+l,c-l) + C'(a+l,6-l,c) + C'(a-l,6,c)- 

Note that the first and the third relations are easily obtained from each other by interchanging the 
first and third component of all dominant weights. Thus 

Xi o X3 and (a, 6, c)o {c,b,a). 

The special recursion relations are obtained from the same products, where some of the com- 
ponents a, b, c of the generic dominant weight take special values 1 and 0. In these cases, the 
decompositions of the three products are modified. To have a complete set of recursion relations, 
every combination of a, b, c with values and 1 needs to be used. 

Solving the following additional recursion relations, we obtain some lowest C-polynomials of A3 
given in Table 9. 

-'^lC'(a,b,l) = 2C(a^b,0) + C'(a+l,b,l) + C(^a-l,b+l,l) + C(a,b-1,2)', 
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XlC^ 



X2C( 


a,6,l) 




a,6,0) 




a,l,c) 




a,l,l) 




1,1,0) 




a,0,c) 




a,0,l) 


X2C( 


a,0,0) 




l,fe,c) 


X2q 


l,b,l) 


X2q 


1,6,0) 


X2q 


l,l,c) 


Xsq 


1,1,1) 


Xaq 


1,1,0) 


Xsq 


l,0,c) 


X2q 


1,0,1) 



a,6,0) = qa+1,6,0) + qa-l,&+l,0) + C(a,b-l,iy^ 

a,l,c) = qo+l,l,c) + qa-l,2,c) + qa,l,c-l) + 2qa,o,c+i); 

a,l,l) = 2qa,l,0) + qa+1,1,1) + C(^a-l,2,l) + ^^(^^,0,2)1 

a,l,0) = qa+1,1,0) + qa-1,2,0) + 2qa,0,l); 

a,0,c) = qa+l,0,c) + C^a-l,l,c) + qa.O.c-l)! 

a,0,l) = 3q„o,0) + qa+1,0,1) + ^(a-1,1,1)' 

a,0,0) = qa+1,0,0) + qa-1,1,0) ; 

i,6,c) = q2,6,c) + 2qo,6+i,c) + qi,6,c-i) + qi,&-i,c+i); 
1,6,1) = 2qi,6,o) + q2,6,i) + 2qo,6+i,i) + qi,6-i,2); 

1.6.0) = 2qo,6+l,0) + q2,6,0) + qi,6-i,i); 

i,i,c) = q2,i,c) + 2qo,2,c) + c'(i,i,c-i) + 2qi,o,c+i); 

1.1.1) = 2qi,i,o) + q2,i,i) + 2qo,2,i) + 2qi,o,2); 
1,1,0) = 2qo,2,o) + q2,i,o) + 2qi,o,i); 

i,o,c) = ^(2,o,c) + 2qo,i,c) + qi,o,c-i); 
1,0,1) = 3^1 + 2qo,i,i) + q2,o,i); 
0,6,c) = qi,6,c) + qo,6,c-l) + qo,6-i,c+i); 
0,6,1) = 2qo,6,0) + + qo,6-i,2); 

0,6,0) = ^"(1,6,0) + qo,6-i,i); 

0,1, c) = 2qo,o,c+i) + qi,i,c) + qo,i,c-i) + qo,o,c+i); 

0,1,1) = 2-'^2 + 3qo,o,2) + qi,i,i); 

o,o,c) = qi,o,c) + qo,o,c-i); 

q2,o,o) + 2^2; ^i-'^2 = qi,i,o) + s-'^s; XiX^ = 4 + qi,o,i); 

2qa+l,6,0) + 2qa-l,6+l,0) + qa,6-l,l) + qa,6+l,l) + qa-1,6,2) + qa+1,6-1,2); 
qa,6-l,0) + qa,6+l,0) + qa-l,6,l) +qa+l,6-l,l); 

2qa,0,c) + 2qa+l,0,c+l) + qa-l,l,c+l) + qa-l,2,c-l) + qa+l,l,c-l) + qa,2,c); 
2qa+l,l,0) + 2qa-l,2,0) + 2qa,0,l) + ^0,2,1) + 2qa+i,o,2) + qa-1,1,2); 
3qa,0,0) + C'(a,2,0) + 2qa+l,0,l) + qa-1,1,1); 
qa,l,c) + qa+l,0,c-l) + qa-l,2,c-l) + qa-l,0,c+l); 
3qa+l,0,0) + 2qa-l,l,0) + C'(a,l,l) + qa-1,0,2); 
qa,l,0) + qa-1,0,1); 

qi,6-i,c) + qi,6+i,c) + q2,6,c-i) + 2qo,6+i,c-i) + 2qo,6,c+i) + q2,6-i,c+i); 
4qo,6+i,o) + 2q2,6,o) + qi,6-i,i) + qi,6+i,i) + 2qo,6,2) + q2,6-i,2); 

qi,6-l,0) + qi,6+l,0) + 2qo,6,l) + q2,6-l,l); 

2qi,o,c) + qi,2,c) + q2,i,c-i) + 2qo,2,c-i) + 2q2,o,c+i) + 2qo,i,c+i); 
4qo,2,o) + 2q2,i,o) + 2qi,o,i) + qi,2,i) + 2q2,o,2) + 2qo,i,2); 
3Xi + qi,2,o) + 2q2,o,i) + 2qo,i,i); 
3qo,o,c+i) + qi,i,c) + q2,o,c-i) + 2qo,i,c-i); 
3q2,o,o) + 4qo,2,o) + 3qo,o,2) + qi,i,i); 
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0,6,c) = C'(o,b_i,c) + C'(o,6+l,c) + ^"(1,6,0-1) + C'(i,b_i,c+1); 
0,6,1) = 2^(1,6,0) + (^(0,6-1,1) + (^(O.b+l,!) + (^(1,6-1,2); 
0,6,0) = C'(0,6+l,0) + C'(l,6-l,l) + C'(0,6-l,0); 
0,l,c) = C'(0,2,c) + C'(l,l,c-1) + 2C(i^o,c+l) + 3C(o,o,c); 
0,1,1) = 2C(i^i^o) + "^^2, + ^"(0,2,1) + 2C(i^2,o); 



0,0,c) 



(0,l,c) 



(i,o,c-i); 





a,6,l) 




a,6,0) 




a,l,c) 


XsC, 


a,l,l) 


XsC, 


0,1,0) 




a,0,c) 




a,0,l) 




a,0,0) 




l,6,c) 




1,6,1) 


X3C 


1,6,0) 


X3C 


l,l,c) 


X3C 


1,1,1) 


X3C 


1,1,0) 


XaC, 


l,0,c) 


X3C 


1,0,1) 


XaC, 


0,6,c) 


XaC, 


0,6,1) 


XaC, 


0,6,0) 


X3C 


0,l,c) 


XzC 


0,1,1) 


XaC( 


0,0,c) 



6 + ^"(0,2,0) + 2C(i^o,i); X2X2, = SXi + C(o,i,i); 

2(^(0,6+1,0) + C(^a-l,b,l) + C'(a+i,6_i,i) + C'(a,6,2); 
C(a-1,6,0) + C'(a+i,b-l,0) + C'(o,b,i); 

2(^(0+1,0,0) + (^(0-1,1,0) + (^(0,2,0-1) + (^(0,1,0+1); 

2C(a,2,0) + 2C(a+l,0,l) + (^(0-1,1,1) + (^(0,1,2); 
3(^(0+1,0,0) + (^(a-1,1,0) + (^(a,l,l); 
C{a-l,0,c) + (^(a,l,c-l) + (^(a,0,c+l); 
2(^(a,l,0) + (^(0-1,0,1) + <^(a,0,2); 
(^(a-1,0,0) + (^(a,0,l); 

(^(l,6,c+l) + (^(l,6+l,c-l) + 2C(o,6,c) + (^(2,6-l,c); 
2(^(1,6+1,0) + 2C(o,6,i) + (^(2,6-1,1) + (^(1,6,2); 
2(^(0,6,0) + (^(2,6-1,0) + (^(1,6,1); 
2(^(0, l,c) + (^(l,2,c-l) + (^(l,l,c+l) + 2(7(2,o,c); 
2(^(1,2,0) + 2(^(2,0,1) + 2C(o,l,l) + (^(1,1,2); 

3(^(2,0,0) + 2^2 + <^(i,i,i); 

2<^(o,o,c) + (^(i,i,c-i); 

2(^(1,1,0) + ^Xs + C(i_o,2); 



(l,6-l,c) + (^(0,6+l,c-l) 



(o,6,c+i); 



2(^(0,6+1,0) + (^(1,6-1,1) + (^(0,1,2); 
(^(1,6-1,0) + (^(0,6,1); 



(l,0,c) 



(0,2,c-l) 



(0,l,c+l) 



(i,o,c); 



(^(0,1,2) + 2C(o,2,0) + 2(^(1,0,1); 



(0,0,c+l) 



+ (^(0,1,0-1); X^ - C(o,o,2) + 2^2- 



E.2 S'-polynomials of A3 

We use the notation S = 'S'(i,i,i) (xi, X2, xa). It is the sum of 24 exponential terms. Generic recursion 
relations are decompositions of the following products, where we assume that a,b,c> 1: 



-^l'S'(a,6,c) = S(^a+l,b,c) + 'S'(a-l,6+l,c) + 'S'(a,6-l,c+l) + S(a,b,c-l)'i 
X2S(^a,b,c) = S(^a,b+l,c) + S (^a+l,b-l,c+l) + 'S'(a-l,6,c+l) + 'S'(a+l,6,c-l) + 'S'(a-l,6+l,c-l) + 'S'(a,6-l,c) ! 
X2,S(^a.,b,c) = 'S'(a,6,c+1) + •S'(a,6+l,c-l) + 'S'(a+l,6-l,c) + S[a-l,b,c)- 



The terms in the decomposition of special recursion relations for S'-functions differ from those of 
C-functions if some <S'(a^ft^c)(^i) ^2; 3:3) = 0. This occurs if some of its dominant weight components 
a, 6, c equal zero. 
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To calculate /S-polynomials explicitly (see Table 9) we use the ^3 characters. The lowest ones 
from the congruence classes # = 0, # = 1,# = 2 and # = 3 are listed below: 



# = 0: 
X(o,o,o) '- 
X(i,o,i) '- 

X(0,2,0) '■ 
X(0,l,2) -- 
X(2,l,0) ■- 

# = 2: 
X(o,i,o) '- 

X{0,0,2) '- 
X(2,0,0) '■ 

X(i,i,i) ■ 



C(o,o,o) = 1, 
3 + ^"(1,0,1)) 

2 + (^(Loa) + C'(o,2,o)) 

3 + 2C(i_o,i) + '^(0,2,0) + C(o,i,2), 
3 + 2C(i,o,i) + (^(0,2,0) + C'(2,i,o); 



(0,1,0) 



^2 



-^2 + C(0,0,2)) 
-^^2 + C'(2,0,0) ) 

4X2 + 2C(o,o,2) + 2C(2,o,o) + C'(i,i,i); 



#=1: 

X(i,o,o) 

X(o,i,i) 

X(2,0,l) 
X(0,0,3) 
X(l,2,0) 

#=3: 

X(OAi) 

X(i,i,o) 

X(l,0,2) 
X(3,0,0) 
X(0,2,l) 



C(i,o,o) = 

2X1 + C(o,-Ll), 

3X1 + C(ojj) + (^(2,0,1)! 

^1 + C(0,1,1) + C'(0,0,3)) 

3X1 + 2C(o,i,i) + C(2fi,i) + C'(i,2,o); 

C(0,0,1) = -^3) 

2X3 + C(i^i^o)> 

3X3 + C(^i_o) + C'(i,o,2)) 

-^3 + C'(1,1,0) + ^"(3,0,0)) 

3X3 + 2C(i^i^o) + ^"(1,0,2) + C'(0,2,l)- 



It should be mentioned that some additional information about polynomials connected with An 
can be found in [39]. We believe that orthogonal polynomials of An are natural n-dimensional 
generalizations of Chebyshev polynomials, but other approaches exist [3]. 



C-polynomials 


# = 


^(1,0,1) 


-4 + X1X3 


C(0,2,0) 


2 — 2X\X'^ + ^2 


C(0,1.2) 


4 - - 2X| + X2XI 


'^'(2,1,0) 


4 - A'lX, + A'2X2 - 2X1 


# = 1 


C(i,o,o) 


Xi 


^(0,1,1) 


-3X1 + 


C(0,0,3) 


3X1 — 3X2 A3 + A3 


C(2,0,l) 


— Ai — 2A2A3 + A^A3 


C(l,2,0) 


5X1 - A2X3 - 2Ai2A3 + AiX| 


# = 2 


C(o,i,o) 


X2 


C(0,0,2) 


-2X2 + xl 


C(2,0,0) 


-2X2 + Xl 




4X2 - 3X1 - 3X1 + X1X2X3 


# = 3 


C(o,o,i) 


X3 


^(1,1,0) 


—3 A3 + A1A2 


C(3,0,0) 


3X3 - 3X1 A2 + Xf 


C(l,0,2) 


— X3 — 2X1X2 + XiX| 


C(0,2,l) 


5X3 — X1X2 — 2X1 x| + X|X3 



^-polynomials 


# = 


5'(2,l,2) 


-I+X1X3 


•^(1,3,1) 


X| — X1X3 


•^(1,2.3) 


1 — xl — X1X3 + X2X3 


'^(:-!,2.i) 


1 + Af A2 - A| - Ai A3 


# = 1 


•^(2,1,1) 


Xl 


•^(1,2,2) 


— Xl + X2X3 


'9(1,1,4) 


Xl — 2X2X3 + X3 


•^(3,1,2) 


— Xl + XfX'3 — X2X3 


•9(2,3,1) 


Xl + X1X2 — X^X3 — X2X3 


# = 2 


"^(1,2,1) 


X2 


•9(1,1,3) 


— X2 + x| 


•9(3,1,1) 


Xl — X2 


•9(2,2,2) 


— X^ — X| + X1X2X3 


# = 3 


•9(1,1,2) 


^3 


•9(2,2,1) 


X1X2 — X3 


•9(4,1,1) 


Xf - X1X2 + X3 


•9(2,1,3) 


— X1X2 — X3 + XiX| 


•9(1,3,2) 


— X1X2 — X3 — XiX| + X|X3 



Table 9: The irreducible C-polynomials and iS-polynomials of A3. 
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E.3 Recursion relations for C- and /S-polynomials of and C3 

The two cases differ in many important respects in spite of the isomorphism of their Weyl groups. 

We write the generic relations for the C-polynomials of the Lie algebras and C3, respectively 
of the simple Lie group 0(7) and Sp{Q). The generic relations for the S'-polynomials are obtained 
by replacing the C symbol by S. 

The variables are denoted by the same symbols Xi, X2, for all algebras of rank 3, namely 
Xj := C(^^ , j = 1, 2, 3. As orbits, they differ for different algebras. They are real- valued for S3 and 
C3, and complex for ^3. 

There are two congruence classes of (01,02,03) for either of the two algebras. 

We have 

^{B^) = 03 mod 2 and #(^3) = 01 + 03 mod 2. 
The congruence numbers add up in a product. For example, in the case of B2, 

#iXlS(^a,b,c)) = #('S'(a,b,c)), #(-^2'S'(a,6,c)) = #('S'(a,&,c))> 

while for C3, we have 

#{^lS{a,b,c)) = #('S'(o,b,c)) + #(^2'S'(a,6,c)) = #('S'(a,6,c))' 

H^{XzS(^a,b,c)) = #iS{a,b,c)) + 1- 

The dimensions of irreducible representations are given by 

dia,b,c){B3) = 7^0(0 + 1)(6 + l)(c + l)(a + 6 + 2)(26 + c + 3)(2o + 26 + 5) X 

(6 + c + 2)(a + 6 + c + 3)(a + 26 + c + 4); 
d{a,b,c)iC3) = 7^(a + 1)(^ + l)(c + l)(a + 6 + 2)(6 + 2c + 3)(a + 6 + 2c + 4) x 

(6 + c + 2)(a + 26 + 2c + 5)(o + 6 + c + 3). 

For 53, we have 

^lC'(a,6,c) = C{a+l,b,c) + C(a-l,b+l,c) + C^a,b-l,c+2) + ^"(0,6+1,0-2) + C^a+l,b-l,c) + C(a-l,b,c)j 

for a, 6 ^ 2, c ^ 3; 

^2C{a,b,c) = C(a,b+l,c) + C(a+l,b-l,c+2) + C'(a-l,fe,c+2) + C'(a+l,&+l,c-2) + (^(0-1,6+2,0-2) + 
C(a+2,6-l,c) + C'(tt+i,b_2,c+2) + C(^a-2,b+l,c) + C'(a_i^(,_i^c+2) + C'(a+i^5_c-2) + 
C(a-l,6+l,c-2) + ^(a^b-l.c)' O ^ 2, 6, C ^ 3; 

^3C(^a,b,c) = C'(a,&,c+1) + C'(a,fe+l,c-l) + C(^a+l,b-l,c+l) + C(a-l,b,c+l) — C^a+l,b,c-l) + 
C(a-l,6+l,c-l) + C(a,b-l,c+l) + C^(a,b,c-l)j a,b,C^ 2. 

For C3, we have the generic recursion relations 

^lC(a,b,c) = C(^a+l,b,c) + C(a-l,b+l,c) + C'(a,b_i,c+1) + C'(a,;,+l,c-l) + 

C(a+i,6-i,c) + ^"((1-1,6,0)) a,b,c^ 2; 

^20(0 ;, ,,) = 5+1 c) + C((j+i c) + C(^ci-l,b,c+l) + C'(a+i,b+l,c-l) + C^a-l,b+2,c-l) + 
C(a+2,b-l,c) + C'(a+i,b_2,c+l) + C'(a_2,6+l,c) + ^(a+i^b^c-l) + C'(a-l,6-l,c+l) + 
C(a-l,6+l,c-l) + C'(a,&-l,c)) O, 6 ^ 3, C ^ 2; 

^3C(^a,b,c) = + C((j^5_,.2,c-1) + C(^a+2,b-2,c+l) + C'(o_2,6,c+l) + ^"(0+2,6,0-1) + 

C(a-2,6+2,c-l) + C'(a,&-2,c+l) + C'(a,&,c-1)) a,b ^ 3, C ^ 2. 
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E.4 C- and S-polynomials of in three real variables 



For group B^, our new coordinates u satisfy the relation = ttj, i = 1, 2, 3. 

The following generic recursion relations for C-polynomials hold true when k,l,'m'^ 2: 

C(k+l,l,m) (^) = '"lC(fe,i,m) {'"')- C(k,l+l,m-2) (^*) -C(it,i-l,m+2) (^*) -C(ifc+l,/-l,m) 

-C(k-l,l+l,m){u)-C(^k-l,l,m){u)-C(^k-l,l+l,m){u)-C(^k-l,l,m)i''J'), 
C(k,l+l,m){u) = U2C(^k,l,m){u)-C!^k+l,l-l,m+2){u)-C(k-l,l-l,m+2){u)-C(k+l,l-2,m^ 

2)(u)-C 

—C(k-l,l+l,m-2) {u)-C(^k+l,l,m-2) {u)-C(^k-2,l+l,m) {u)-C(^k+2,l-l,m) (^) > 
C(k,l,m+1) (^) = '"3C(fe,i,m) (^) -C'(fc+l,/-l,m+l) (^) -C'(ifc,/-l,m+l) (^) -(^(jfc-l./.m+l) (■") 
-C(k-l,l+l,m~l){u)-C(^k,l+l,m-l){u)-C(^k+l,l,m-l){u)-C(^k,l,m-l){u) 

Remaining recurrence relations except for the lowest polynomials are listed below: 

C{k+i,i,o){u) = uiC(^k,i,o)iu)-C(^k-i,i+i,o)iu)-C(^k+i,i-i,o){u)-C(^k-i,i,o){u)-C(^k,i-i,2)iu), 
C(k,i+i,o){u) = U2C(^k,i,o)iu)-C(k-2,i+i,o)iu)-C{k+2,i-i,o){u)-Ci^k-i,i,2){u)-C(^k,i-i,o)iu) 

-C'(fc+l,«-l,2)('")-C'(fc-l,«-l,2)(^)-C'(fc+l,/-2,2) (u), 

C{k+i,o,o){u) = uiC(^k,o,o)iu)-C(^k-i,i,o)iu)-C{k-i,o,o)iu), 

C(0,/+1,0)N = «2C(o,Z,o)(«)-C(2,«-l,0)(«)-C(0,Z-l,0)N-C(l,«-l,2)W-C(l,Z-2,2)N, 
C(0,0,m+1)(^) = ■"3C'(0,0,m)('")-C'(0,l,m-l)('")-C'(l,0,m-l)('")-C(0,0,m-l)(^^)) 
C(k+l,0,m){u) = UlC!^k,0,m){u)-C(^k,l,m-2)iu)-C(k-l,l,m){u), 

C(k,0,m+l)(,u) = U3C(^k,0,m)iu)-C(^k-l,0,m+l){u)-C(^k-l,l,m-l)iu)-C(^k,l,m-l)(,u) 

~C(k,Q,m-l) ('")~C'(fe+l,0,m-l)(^)) 
C(0,;+l,m)('") = ■"2C(o,/,T„)('")-C(i^;_i^^+2)(^)-C(l,/-2,m+2)('")-C(0,«-l,m)('") 

-C(2,/-l,m)(^*)-C(l,/+l,m-2)(^t)-C(l,i,m-2) 
C(o,l,m+l){u) = 1t3C(o,i,m)('")-C(l,i-l,m+l)(^*)-C(0,/-l,m+l)('")-C(0,;+l,m-l)('") 

-C'(0,i,m-1)(^)-C'(l,i,m-1)(^)- 



The lowest C-polynomials of B^ were calculated explicitly and arranged in Table 10. 
As in the previous case, we can use the Weyl character formula or generic recurrence relations 
for ^'-polynomials of B^ valid for k,l,m ^ 2: 

^(k+l,l,m) (u) = UiS(^k,l,m) (^) -^(k,l+l,m-2) {u) -S!^k,l-l,m+2) (u) -S(^k+l,l-l,m) (^) 
-'S{k-l,l+l,m) {u) -S(^k-l,l,m) {u) -S(^k-l,l+l,m) {u) -S(^k-l,l,m) {u) , 

'S(k,l+l,m) (u) = U2S(^k,l,m) {u) -'S(^k+l,l-l,m+2) {u) -S(^k-l,l-l,m+2) {u) -S(^k+l,l-2,m+2) {u) 
—^{k+l,l,m-2) {u) ->5(fe_2,;+l,m) (^) -^{k-l,l,m+2) {u) -S(^k-l,l+2,m-2) {u) 
~^{k+l,l+l,m-2) —^{k-l,l+l,m-2) ("") —^{k+2,l-l,m) (^) > 

^(k,l,m+l) ("") = '"3'5(fe,i,m) (u) -S(^k+l,l-l,m+l) (^) -^(k,l-l,m+l) (^) -^{k-l,l,m+l) (^) 

-'S(k-l,l+l,m-l) iu)-S(^k,l+l,m-l) {u) -S(^k+l,l,m-l) {u) -S(^k,l,m-1) (^) • 
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Table 10: Lowest C-polynomials of B3 split into two congruence classes # = and # = 1. 
E.5 C- and /S-polynomials of C3 in three real variables 

Low order C-polynomials of group C3 are listed in Table 11. Higher-order polynomials can 
obtained from the recurrence relations. Generic recursions for k,l,m 2 are: 

C(k+l,l,m){u) = UiC(^k,l,m)iu)-C{k,l-l,m+l)iu)-C(k+l,l-l,m)iu)-C(^k-l,l+l,m){u) 

— C(k-l,l,m){'u) — Ci^k,l+l,m-l)(fJ-)i 

C{k,l+l,m)iu) = U2C(k,l,m)iu)-Ci^k+l,l,m-l)iu)-C\k-l^l^rn+l)iu)-C\k^l-l^m)iu) 

~C'(A:+l,/+l,m-l) (''-')~C(fc-2.i+l,m) — C(fc-lJ-l,m+l)('") — C(fc-lJ+l,m-l) 

— (^(A:+2,/-l,m)('") — C^(fc+l,/-l,m+l)('") — C'(fc+l,i-2,m+l)('") — C'(A:-l,«+2,m-l)(^^)) 
C{k,l,m+1) (u) = U3Ci^k,l,m) iu)-C\k,l^2,m+l) iu)-C\k^2,l,m+l) ("") -C'(fe,«+2,m-l) (^) 

—C(k,l,m-l) {u)-C(^k+2,l-2,m+l) {u)-C(^k-2,l+2,m-l) {u)-C(^k+2,l,m-l) {u) ■ 

Additional recursions: 

C(fc+i,o,o)(^) = ^iC(fc,o,o)M-C(fc-i,i,o)(^^)-C(fc-i,o,o)(^)' k>l; 
C(k+i,i,o){u) = uiC(^k,i,o){u)-C(^k,i-i,i){u)-C(^k-i,i,o)(.u) 

-C(k+i,i-i,o)i'<J')-C(k-i,i+ifi){u), k,l > 1; 

C(k,l+lfl)iu) = U2C(^k,l,0)iu)-C\k-l,l,l)iu)-C(^k,l-l,0){u)-C(^k+2,l-l,0)iu)-C(^k-2,l+l,0)i'i^) 

-C^k-i,l-i,i)iu)-C^k+i,i-i,i)iu)-C(^k+i,i-2,i){u), k,l>2; 

C(0,Z+l,m)(«) = ''^2C'(o,/,rn)('")-C'(l,i,m-l)('")-C'(0,i-l,m)(^)-Q2,Z-l,m)(«)-C'(l,Z+l,m-l)('") 
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-C(l,/-l,m+l)H-C(l,/-2,m+l)(^)> ^ > 2, m > 1; 
C(0,;,m+1)('") = 1*3C(0,i,m)('")-C(0,/-2,m+l)(^)-C(0,i,m-l)(^*)-C(2,;-2,m+l)('") 

-C(2,/,m-l)W-C(0,i+2,m-l)W> ^ > 2, m > 1; 
C(fe+l,0,m)(^*) = UlC(^kfi^jn){u)-C(^k-l,l,m){u)-C(^k,l,m-l)iu)-C(^k-l,0,m){u), k,m > 1; 
C(fe,0,m+1)(^) = '"3C(jk,0,m)('")-C'(fe-2,0,m+l)(^)-C'(it,0,'n-l)('")~^{fc-2,2,m-l)(^) 

-C(ifc+2,0,m-l)('")-C(ifc,2,m-l)(^)> ^ > 2, m > 1; 

C(o,Hi,o)W = '"2C(o,i,o)('")-C(o,i-i,o)W-C(2,i-i,o)W-C(i,/-i,i)('")-C(i,i-2,i)W, ^ > 2; 

C(0,0,m+1)('") = ■?^3C(0,0,m)(^*)-C(0,0,m-l)(^)-C(2,0,m-l)(^*)-C(0,2,m-l)(^)) m > 1. 

Generic recursions for S'-polynomials hold true when k,l,m^2 

-'S{k-l,l,m){u)-<S{k,l+l,m-l)iu), 
^(k,l+l,m) (■") = U2S(^k,l,m) («) -'S{k+l,l,m-l) i'^) -^{k-l,l,m+l) {u) -«5(fc+2,i-l,m) (^*) 

—^{k-2,l+l,m){u)-S(^k_l^l_l^rn+l){u)-S(^k-l,l+l,m-l) {u)-S(^k_^_i^i_^_i^rn-l)iu) 

S{k+l,l-l,m+l) {u)—S(^k+l,l-2,m+l) >5(fc_i /+2,m-l) {'"')~^ik,l-l,m) (^)) 

^(k,l,m+l) (■") = U3S(^k,l,m) (u) S(k,l-2,m+l) {'"') -<S(k-2,l,m+l) (u) -<S(k,l+2,ni-l) ("") 

—^{k+2,l-2,m+l) (■") —^{k-2,l+2,m-l) (■") —^{k+2,l,m-l) {u) -S(^k,l,m-1) (^) • 
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